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We reinterpret the dilepton emission rate from a hadronic gas expected to be 
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produced in heavy ion collisions in terms of the spectral functions available 



4^3 ' from the e + e annihilation and the r lepton decays experiments. We take 

Qh into account all possible hadronic state especially the multi-pion contributions 



to the dilepton emission and the parity mixing phenomenon due to the soft 



final-state corrections. A new compilation of the experimental data for the 

h : 

spectral functions is presented. 
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The lepton pair emissions are the ideal probes of the dense system of strongly interacting 
particles expected to be produced in ultrarelativistic heavy nucleus-nucleus collisions [|TJ. 
Once produced, they will decouple from the strongly interacting system without further 
interactions and thus carry the information on the dynamical properties of the hot nuclear 
system. The theoretical predictions on the emission rate, both from a quark-gluon plasma at 
high temperature (T) and from a hot hadron gas at relatively lower T, have been extensively 
discussed in literatures [@-|6|]. There already exist some experimental data from various 
CERN SPS experiments (for a recent review, see @), and more thorough investigations will 
be made at upcoming RHIC and LHC experiments. The quantitative interpretation of the 
data is beginning to be possible [§,[| 

The main theoretical uncertainties in calculating the dilepton emissions rate come from 
the fact (apart from the crucial assumption that the system thermalizes in short time and 
stays in local equilibrium) that the system is strongly interacting and the perturbation 
theory fails. In a hot hadronic gas, though the effective degrees of freedom are mainly pions, 
pions interact strongly among themselves in many resonance channels, and effectively the 
resonances have to be taken into account. One obvious example is the pion electromagnetic 
form factor which exhibits the p-meson dominance. Near the resonances the naive pole 
approximations often violate the unitarity of the partial waves. There may also be some 
potential double counting problems associated with the emission rate at the location of the 
resonance peaks, as to whether or not the real resonance dilepton decays should be counted 
separately. Also, the multi-pion state emissions, e.g. 3tc — > £ + £~ and Air — > £ + £r may not be 
so suppressed as they can proceed through the resonance poles, e.g. 3n — > u*, <fi* — > £ + £~ and 
47r — > 7r<2i — > £ + £~ . Taking all these reactions into account by the kinetic theory is a difficult 
task though not impossible. One may always wonder the validity of a calculation and if 
the omission of certain reactions oversimplifies the life. More importantly, process involving 
hadrons in the final state such as nir — ► £ + £~mTT may have been entirely overlooked. Even 
in the phase of quark-gluon plasma where quarks and gluons are supposed to be weakly 
interacting, the high order corrections in a s may not be so small. 
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The aim of this paper is to use and to develop the expression that allows the non- 
perturbative treatment of the strong interactions and avoids the detailed counting of the 
relevant reactions. Such an expression is known to exist Pfl|T0|| in the form of a low temper- 
ature expansion when T 2 is smaller compared to the dilepton invariant mass squared M 2 . 
We shall reinterpret the dilepton emission rate in terms of the spectral functions of some 
hadronic currents that are known experimentally in the low energy e + e~ annihilation process 
and t lepton decays. While we shall comment on its application to the quark-gluon plasma, 
our approach is most relevant for the case of a hadronic gas. The recent development in 
calculating the correlation functions of vector currents in a low temperature pion gas using 



the PCAC and current algebra |TJ] enables us to obtain the important corrections to the 
rate with soft pions emitted together with the dilepton in the final state. Our parameteri- 
zation clearly exhibits the mixing phenomenon of vector and axial vector currents at finite 
temperature and its relation to the dilepton production. We present a new compilation of 
the experimental data for the spectral functions in the energy range from 360 MeV to 2150 
MeV, relevant to the dilepton yield from a hadron gas. 

In the lowest order in the electromagnetic interactions, the emission rate can be expressed 
non-perturbatively in terms of the correlation functions of the hadronic electromagnetic 
current. When the lepton mass is ignored, the rate reads [p|,|3|.|l0f 

where 

W^(q) = I d*xe-^((J™(x)Jt™m . (2) 



q is the time-like 4-momentum of the lepton pair, J^ m is the electromagnetic current for 
hadrons and ((•■■)) stands for the thermal average at a give temperature T. Note that 
W flv (q) is the thermal structure function, and there is no time ordering in the product of 
J^ m (x) and Jt em (0), which is crucial in relating W^ u directly to the spectral functions, flip 
is shown to agree with the relativistic kinetic theory reaction by reaction if one inserts the 
complete set J2f = 1 f° r each component \I) of the thermal density matrix 



= EE / d*e-^(I\J™(x)\F)(F\4™(0)\I)^T- > ( 3 ) 

F I ^ 

where Z is the canonical partition function of the system. (^[) describes the all possible 

transitions I — > F£ + £~. Summing up the disconnected parts of the matrix elements, one 

recovers the appropriate Fermi or Bose distribution functions f{Ej) = [eP Ei ± for each 

initial-state particle and the suppression or enhancement factor (l=p/(£ , i )) for each final-state 

particle. W^q) has no non-trivial zero temperature limit. To get the leading temperature 

dependence, one alternatively sums up all initial states |/)'s using the completeness Y^i \I){I\ 

and the energy conservation Ej = Ep + q° [[TIJ 

r p-PE F 
W ltv (q)=e-M > Y I J d A xe^(F\Jl™(x)J™(0)\F)^- 

= J d 4 xe^((jt em (x)J-(0))) = esprit) ■ (4) 

p^v(q) now has a non-trivial limit as T — > 0. Though \F)'s and |J)'s seem to switch their 
positions in ([|) and the same kinetic theory result would be obtained if one uses the 
identity (1 =F f{E))e~ l3E = f(E). 

We would like to compute p^ u (q) for a low temperature hadronic gas, not reaction by 
reaction as one does in the kinetic theory, but to obtain a low T expansion, in which 
each term is the summation of all possible reactions. The leading term is independent of 
temperature, which is obtained by taking the limit T —>■ where only the vacuum state in 
|F)'s contributes to the connected part, namely 

PM = / d 4 xe^(0\ Jf*(x) J7(0)|0) . (5) 

When all possible initial states |J)'s are saturated between the two currents, (|5|) corresponds 
to the "hard" particle contributions illustrated in || where "hard" means "not necessarily 
soft". As we are interested in the kinematical region T 2 <C q 2 , any non- vacuum states 
|F)'s are necessarily soft since they are kinematically limited to E < T by the presence of a 
thermal distribution function. The low but finite temperature corrections to the leading term 
p e J£(q) arise from the soft final-state particle emissions and the soft initial-state particles, 
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both are limited to E < T. In a hadronic gas, the soft modes are the pions. One may use 
the soft pion theorem based on the PCAC to calculate the soft pion corrections. Unlike the 
case in the quark-gluon plasma where there is a contribution of order T 2 , the insertion of 
the soft pions (the initial states) in the intermediate states of the correlator does not pick 
up a contribution of order T 2 . In fact, the extra powers of soft thermal pion momentum 
leads to a suppression factor of order T 4 /M 4 in the chiral limit, as shown by Leutwyler and 
Smilga in [12|] and Eletsky in fil l. This is the consequence of the low temperature theorem 



|T2|1 (the finite temperature version of the Adler theorem) that the chiral symmetry protects 
all of the masses from picking up a contribution of order T 2 . We shall ignore the effects of 
the initial-state soft corrections with the understanding that they start at 0(T 4 /M 4 ). 

On the other hand, the corrections from the soft final-state thermal pions, as shown by 



Dey, Eletsky and Ioffe leads to the mixing of vector and axial vector correlators. The 
essence is to approximate the matrix elements containing soft pions with their threshold 
value at zero pion momentum, which is completely determined by the current algebra or the 



symmetry of the system |T3[ . For example 

1 

F 2 



mr+l J] cm (x) j; m (0)\nn + ) = ~({n - l)n + \J^(x)V 3 (0) - J^{x)A+{0) 



-A+(x)A-(0) + V*J™(0)\(n - l)7r+> (6) 
= ^(-irS^OlK 3 W 3 (0) - Al(x)Al(0)\0) , 

where V 3 and A? are the third (neutral) component of the vector- and axial- vector- isovector 
currents, F n = 93 MeV. By doing so, we may compute all reactions with soft pions in the 
final state / — > £ + £~nir . It is worth pointing out that although J^ m contains all flavors 
(u, d, s, c, b, t), the soft pion emission selects only two flavors (u, d) since and A 3 ^ contains 
no strange or higher flavor contents. Also, when only u, d flavors are effective at lower energy 
range, J^ m = V 3 + where is the baryonic current which contains both isovector 

(/ = 1) and isosinglet (/ = 0) parts, while V 3 and A 3 ^ are purely isovector. When n > 2, 
(^) combining with the integration of the Bose distribution functions gives the corrections 
of order (T 2 /6F 2 ) 2 or higher. However, the multi-pions in the final state can interact with 



themselves and cause further phase shifts. Inclusion of the tttt rescattering effects gives a 
correct T 4 term |IT|. We neglect the non-zero pion momentum that causes a correction 



suppressed by M 2 and summarize the corrections to T 4 in the spectral tensor 

PM = PTM - (c - j)(pJM - ffc{q)) , (7) 

where e = T 2 /6F 2 . p^ u (q) and p^ u (q) are the correlators defined in (|) with J° m replaced 
with V 3 and A 3 ^ respectively. Clearly, the dilepton production from a thermal source comes 
not only from the vector channel but also from the axial vector channel since they can 
interchange their parity identities through the interaction with a thermal background. The 
soft corrections tend to suppress the contributions from the vector channel and to increase 
the rate by adding to it the axial channel contributions. 

Since J° m and V? are conserved currents (neglecting the isospin breaking), p e ^(q) and 
Pta/io) have only transverse (spin one) part. receives both spin one and spin zero (pions) 
contributions. However, the longitudinal part (oc q^q y ) does not contribute when contracted 
with the lepton pair tensor. Define the spectral functions p^'(s) (i = em, V, A) 



p«(g) = 2v{q,q u - q 2 g, v )9(q°)P { Hs) , (8) 
where s = q 2 , one obtains the dilepton emission rate in terms of three spectral functions 

^ = ^-^{^ m W-(^-y)lA' v W-A«)]} , 0) 

or if the transverse mass and the rapidity are integrated over 

dR 4a 2 , f _ , , _ 4 , e 2 ^ 



(np ^-MTK X {M/T) |p cm (M) - (c - -)[p v (M) - p A (M)]| . (10) 

As we have noted earlier, because the spectral tensor does not involve the time ordering of 
the currents, it is directly represented by the spectral functions on the physical axis q 2 > 0. 
In the case of time ordering correlator, the spectral function is (twice) the imaginary part 
of the correlator. To get the real part, one has to perform the dispersion integral of the 
spectral function. Clearly, the usual pole-plus-continuum approximation widely used in the 
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Weinberg's sum rules QT3j] would not do since it is precisely the energy dependence of the 
spectral densities that determine the shape of the dilepton spectrum. 

Other than a theoretical estimate, a more sound procedure is to use experimental date 
which have been accumulated over years and measured accurately enough to truly represent 
these spectral functions. Obviously, p em (s) can be extracted from the total hadronic cross 
section (the R value) of e + e~ annihilations (note that <r(had — > e + e~) = 4<r(e + e~ — > had) 
because of the helicity summation) 

n em/ c x _ sa(e + e- -> had) _ R(s) 

9 K ) ~ \^a 2 - i 27r 2 • W 

To get p v (s), we have to first restrict the open flavors to u, d. This can be done effectively by 
selecting the events in which only pions are produced. We also need to eliminate the 1 = 
events in e + e~ annihilations. Since the 1 = states will have the G parity G = C(-l) 1 = —1, 
which can only decay into an odd number of pions (a pion has G = —1). So by selecting 



the events consisting only of even number of pions, we have for p [0.0 

p A can be extracted from the differential probabilities of the r lepton decays into odd number 
of pions r — > u T (2n + l)n [ 16| 



A/ \ S7im l v dT(r -> u T (2n+ 1)tt) 

P{ ) G 2 F cos 2 e c (m 2 T + 2s)(m 2 T -s) 2 ^ ds ' K ' 

where s is understood to be the invariant mass squared of (2n + l)7r's, Gp is the Fermi 
constant and 8 C is the Cabibbo angle. 

We present a new compilation of the data for those spectral functions in the energy 
range fro 360 MeV and 2150 MeV where the data are rich. The total hadronic cross section 
of e + e~ annihilation above \/s = 1420 MeV is measured by Adone 772 and MEA |T 



below which the ORSAY data exist [18] but are fairly sparse, especially at p,uj,cj) peaks. 



We fill in more data from OLYA and CMD e + e — > 2tt measurement |TjJ, Orsay DM1 



e + e — > 3tt measurement [20|, and VEPP-2M e + e — > K + K ,K$Kl |2^] in the resonance 



region. To extract p , we consider, in addition to e + e~ — > 27r, the multi-pion reactions 
e + e~ — ► 27T + 27r _ , 7r + 7r~27r , 27r + 27r~27r°, 37r + 37r~ which are available from Adone P2"| . For 
p A , we assume the equality of T(r — > z/ r 7r + 27r _ ) and T(r — > z/ T 7r _ 27r°) which is expected if 
the 37r state is dominated by the ai(1230) resonance, and use the recent ARGUS data ||23|| . 
Clearly, p A thus obtained is limited to y/s < m T ~ 1.7 GeV. However, the rapid fall off of 
p beyond m ai indicates a rather negligible contribution at a larger yfs. 

We plot our compilation of the numerical values of the three spectral functions in Fig. 
1 for the convenience of further phenomenological studies. Data with poor statistics are 
removed by hand. The dilepton emission rate at a typical temperature T = 160 MeV is 
plotted in Fig. 2 for the invariant mass ranging from 360 MeV to 2150 MeV. The dotted line 
is the leading contribution without the soft final-state emissions. An enhancement located 
around m ai ~ 1230 MeV is visible due to p A . 

We have reinterpreted the dilepton emission rate in terms of spectral functions available 
directly from the data. The leading contribution is related to the total inclusive hadronic 
cross section of e + e~ annihilation, which sums up all possible combinations of mesons and 
resonances such as P + P, V + V, V + P — > £ + £~ recently considered by Gale and Lichard [[5], 
and three- and four- body combinations, and many more. It is free of theoretical uncertain- 
ties in calculating these non-perturbative processes and avoids the possible double counting 
of the light mesons and their resonances. As emphasized by McLerran and Toimela |§, the 
leading term by no means contains only the hard particle contributions. It is merely not 
necessarily soft which includes the soft particle contributions as well. The parity mixing 
phenomenon, also considered in PJ9|,|24|| is explicit in our reinterpretation by the presence of 
the axial spectral density p A , which arises from the soft final-state emission of pions. Our 
formalism is applicable to a hadronic gas as long as T 2 <C M 2 , T 2 -C 6F 2 and T < T c . It 
may be possible to extend the approach to the case of quark-gluon plasma for the high in- 
variant mass region M > 2 GeV, where a(e + e~ — > had) ~ a(e + e~ — > qq) and a s corrections 
to all orders would be automatically included . When M < 2 GeV, the parton degree of 
freedom are not apparent in a(e + e~ —* had). The existence of the parton degrees freedom 
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at T > T c even for low invariant mass parton pair is a fundamental assumption. To calculate 
the actual dilepton multiplicity in order to comapare it with the experimental data, one has 
to integrate the rate over the space-time history of the thermal source, which involes more 
assumptions and approximations, and will be reported elsewhere. 

I would like to thank Y. Kluger, A. Kovner and X.-N. Wang for useful discussions. This 
work was supported in part by the Director, Office of Energy Research, Office of High Energy 
and Nuclear Physics, Division of High Energy Physics of the U.S. Department of Energy 
under Contract DE-AC03-76SF00098 and by the Natural Sciences and Engineering Research 
Council of Canada. 
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Figure Captions 



Fig. 1 Spectral functions p cm (s), p V (s), and p A (s) from e + e annihilation and r decays. 
Note that p cm (s) extends beyond the vertical scale at the ui and resonance peaks. 

Fig. 2 Dilepton emission rate from a T = 160 MeV thermal source. Dotted line (error 
bars not shown) represents the leading temperature contribution without the soft pion 
emissions. The error bars are calculated assuming the statistical independence of 
different experiment sets. 
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